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Abstract 


This Letter presents a search for a heavy neutral particle decaying into an opposite-sign 
different-flavor dilepton pair, or using 20.3 fb”* of pp collision data at 

= 8 TeV collected by the ATLAS detector at the LHC. The numbers of observed can¬ 
didate events are compatible with the Standard Model expectations. Limits are set on the 
cross section of new phenomena in two scenarios: the production of Vt in 7?-parity-violating 
supersymmetric models and the production of a lepton-flavor-violating Z' vector boson. 
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Lepton-flavor-violation (LFV) in the charged sector, if observed at present sensitivities, would be a clear 
signature of new physics. Such signatures occur in several new physics scenarios, including /?-parity- 
violating (RPV) [1] supersymmetry (SUSY) [2-10] and models with an additional heavy neutral gauge 
boson Z' [11] allowing for LFV couplings. 

In RPV SUSY, the lagrangian terms allowing LFV can be expressed as ^AijkLjLjek + A'.jj^LiQjdk [1], 
where L and Q are the S U{2) doublet superfields of leptons and quarks; e and d are the S U{2) singlet 
superfields of leptons and down-like quarks; A and A' are Yukawa couplings; and the indices i, j and 
k denote fermion generations. A r sneutrino (y,-) may be produced in pp collisions by dd annihilation 
and subsequently decay to ep, er, or pr. Although only Vr is considered here in order to compare with 
previous searches performed at the Tevatron, the results of our analysis apply to any sneutrino flavor. 

The Sequential Standard Model (SSM), where the Z' boson is often assumed to have the same quark and 
lepton couplings as the Standard Model (SM) Z boson, can be extended to include LFV couplings for the 
Z'. The Z' ep, er, or pr couplings (Qu, Qn, or Q 23 ) [12] are typically expressed as fractions of the 
SSM Z' {€ = e,p, t) coupling. 

The CDF [13], DO [14], and ATLAS [15] collaborations have searched for a Vr in LFV final states and 
placed limits for various Vr mass hypotheses. Both the CDF [16] and ATLAS [17] collaborations have 
placed limits on Q \2 as a function of the Z' mass. 

This Letter describes a search for a neutral heavy particle (vt or Z') decaying into e-p'^ {ep), e~T'^^^ (er), 
or p~T^^^ ipT) using pp collision data collected at = 8 TeV, where Thad is a r lepton that decays into 
hadrons. The ATLAS detector is described in detail elsewhere [18]. Events are selected with a three-level 
trigger system that requires one or two leptons {e or p) with high transverse momentum (pt)- The dataset 
has a total integrated luminosity of 20.3 ± 0.6 fb“^ where the uncertainty is derived following the same 
methodology as that detailed in Ref. [19]. 

Electrons are required to have pj > 25 GeV, \p\ < 1.37 or 1.52 < Ipl < 2.47 [20], and satisfy the 
“tight” selection in Ref. [21], which was modified in 2012 to reduce the impact of additional inelastic 
pp interactions, termed pileup. Muon candidates must have py > 25 GeV, Ipl < 2.4 and be reconstructed 
in both the inner tracker detector and the muon spectrometer. The muon momenta measured by the inner 
detector and muon spectrometer must match within five standard deviations of their combined uncertainty. 
Good quality reconstruction and pT resolution at high momentum are ensured by requiring a minimum 
number of associated hits on the inner detector track [22] and in each of the three muon spectrometer 
stations. 

Candidate events must contain at least one primary interaction vertex reconstructed with more than three 
associated tracks with pj > 400 MeV. If there is more than one such vertex, the one with the highest 
sum of p| of associated tracks is chosen. The longitudinal impact parameter is required to be smaller 
than 2 mm for candidate electrons and smaller than 1 mm for candidate muons. It is further required 
that the transverse impact parameter is less than six times its resolution for candidate electrons, and that 
the transverse impact parameter is smaller than 0.2 mm for candidate muons. A calorimeter isolation 
criterion lEj < 0.06 and a tracker isolation criterion p^^^° ^/pT < 0.06 are applied for both the 

electrons and muons, where is the transverse energy deposited in the calorimeter within a cone 

of size AT? = 3 j{Ar])^ + (A^)^ = 0.2 around the lepton, and is the sum of the px of tracks with 

Pt > 1 GeV within a cone size of 0.4 around the lepton. Ej and px are the lepton transverse energy and 
momentum, respectively. 
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Hadronic decays of r leptons are characterized by one or three charged tracks associated to a narrow 
energy cluster in the calorimeters [23]. This search uses Thad candidates with only one charged track due 
to reduced identification and reconstruction efficiency for high-px t decays wifh fhree charged hacks. 
Boosfed-decision-hee mulfivariafe discriminators, based on hacking and calorimefer informafion, are 
used fo rejecf jefs and elecfrons misidenfified as Thad- The Thad candidafes musf have \r]\ < 2.47 and 
> 25 GeV. Candidafes wifh \r]\ < 0.03 are removed fo exclude a region wifh increased misidenlificalion 
of elecfrons due fo reduced coverage by fhe inner defecfor and calorimefers. 

Jefs are reconsfrucfed from clusfers of energy in fhe calorimefer using fhe anfi-k( algorifhm [24] wifh 
radius paramefer R - 0.4. Jef energies are calibrafed using Monfe Carlo (MC) simulafion and a combina- 
fion of several in sifu calibrafions [25]. The calculafion of fhe missing hansverse momenfum vecfor p 
(wifh magnifude S™®®) is based on fhe vecfor sum of fhe calibrafed px of reconsfrucfed jefs, elecfrons, 
and muons, as well as calorimefer energy clusfers nof associafed wifh reconsfrucfed objecfs [26]. 

Candidafe signal evenfs are required to have exacfly fwo leptons, of opposite charge and of differenl 
flavor, salisfying fhe above lepfon selection criteria. The fwo lepfons are required to be back-fo-back in 
fhe azimufhal plane wifh > 2.7, where is fhe cf) difference befween fhe fwo lepfons. Due fo 

fhe presence of fhe undefecfed neufrino in fhe t decay, fhe ^x of fhe Thad candidafe is required fo be less 
fhan fhe ^x (px) of fhe elechon (muon) for fhe CThad (p^had) channel. 

A collinear neufrino approximafion is used fo reconsfrucf fhe dilepfon invarianf mass {mu>) in fhe CThad 
and pThad channels. In fhe hadronic decay of a t lepton from a heavy resonance, fhe neufrino and fhe 
resulfanf jef are nearly collinear. The four-vecfor of fhe neufrino is reconsfrucfed from fhe p™®* and p of 
fhe Thad jot- Four-vectors of fhe elecfron or muon, Thad candidafe and neufrino are fhen used fo calculate 
mu’. For CThad and pThad signal evenfs, fhe above fechnique significanfly improves fhe mass resolution 
and search sensifivify. 

Evenfs wifh niw< 200 GeV form a validafion region fo verify fhe background modeling, and evenfs wifh 
mw > 200 GeV are used as fhe search region. 

The SM processes fhaf produce final sfafes can be divided info fwo cafegories: processes fhaf pro¬ 
duce fwo prompf lepfons such as Z/y* —> tt, tt, single-fop Wt channel, diboson producfion, and processes 
where one or more photons or jefs are misidenfified as lepfons, predominanfly IT/Z + y, IT/Z-i-jefs, and 
mulfijef evenfs. The decay of a t fo an elecfron or a muon is considered as prompf producfion. For fhe 
^Thad (p^had) channel, additional background can originafe from fhe Z/y* ee (pp) process if one lepfon 
is misidenfified as a Thad candidafe. The confribufions of fhese processes are even larger wifh respecf fo 
fhe Z/y* —> TT background, since fhe final sfafes e or p are usually harder fhan fhose from lepfonic t 
decay. 

Confribufions from processes in fhe prompf fwo-lepfon cafegory, as well as phofon-relafed and Z/y* —» 
ee ijufj.) backgrounds, are esfimafed using MC simulafion [27]. The defecfor response model is based 
on fhe geant4 program [28]. Lepfon reconsfrucfion and idenfificafion efficiencies, energy scales, and 
resolufions in fhe MC simulation are correcfed fo fhe corresponding values measured in fhe dafa. Pileup 
is included to mafch disfribufions observed in fhe dafa. Top quark producfion is generated wifh mc@nlo 
v4.06 [29] for tt and single-fop, fhe Drell-Yan process {Z/y* —> is generated wifh alpgen v2.14 [30], 
and diboson processes are generated wifh herwig v6. 520.2 [31]. Samples of Wy and Zy evenfs are 
generafed wifh sherpa v1.04 [32]. These generafed samples are normalized fo fhe mosf accurafe available 
cross-section calculafions. For fhe dominanf backgrounds, fhe Drell-Yan processes are correcfed fo nexf- 
fo-nexf-fo-leading order (NNLO) [33], and tt is correcfed fo NNLO, including sofl-gluon resummation to 
next-to-next-to-leading-logarithm order [34]. 
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Since it is difficult to model misidentilication of jets as leptons, particularly at high pj, the VK+jets and 
multijet backgrounds are determined from control regions in the data. The IT+jets background is determ¬ 
ined in a control region selected with the same criteria as used for the signal selection except requiring 
> 30 GeV (to enhance the W contribution) and requiring that the electron or muon pj be less than 
150 GeV (to eliminate potential signal). Simulation studies indicate that there is negligible multijet back¬ 
ground in this control region. For the CThad and prhad channels, the number of events in the control region 
is corrected for the other SM background sources using MC samples. For the ep channel, the number of 
VF-i-jets events in the control region is too small to yield a statistically meaningful measurement. Instead, 
the control region is enlarged by removing the isolation criterion on one lepton, and the VF-i-jets contribu¬ 
tion is estimated using the lepton distribution to fit the data with the MC predictions for other 

SM processes (dominant at low values of the isolation variable) and IF-i-jets (dominant at large values). 
For the IT-i-jets background in all three channels, the extrapolation factor from the control region to the 
signal region and the shape of the mtp distribution are taken from the IT-i-jets MC sample. 

The contribution from multijet production is estimated from a control region with the same selection as the 
signal region except that the leptons are required to have the same electric charge, under the assumption 
that the probability of misidentifying a jet as a lepton is independent of the charge. The number of events 
in this same-sign control region is corrected for contributions from backgrounds with prompt leptons and 
from VF-rjets backgrounds using the procedure described above. The assumption of charge independence 
of the jet misidentification rates is tested in a multijet-enriched region with two nonisolated leptons. 
After subtraction of other SM backgrounds, the ratio of opposite-sign to same-sign events is found to be 
1.07 ± 0.06 (stat.) ± 0.02 (syst.). 

The background estimates are verified in validation regions defined with the same selection criteria as 
for the signal regions but with 1.0 < |A0/';/| < 2.7. In these validation regions, simulation studies show 
the backgrounds have compositions similar to the signal regions, and the predictions agree with the data 
within 20% over the entire mw range. An uncertainty of 20% is hence placed on the total background 
estimate that is used in the final results. 

MC signal events are generated with herwig v6.520.2 for and pythia v8.165 [35] for Z'. Samples 
are produced with Vt and Z' masses ranging from 0.5 to 3 TeV. Signal cross sections are calculated to 
next-to-leading order for y^, and NNLO for Z'. The theoretical uncertainties are taken from an envelope 
of cross-section predictions using different parton distribution function (PDF) sets and factorization and 
renormalization scales [36,37]. 

The selection efficiency, including t decay branching ratio if r is involved, for niy^ - 2 TeV are 42%, 
14%, and 10% in the ep, ct and pT channels, respectively. The corresponding numbers for a Z' boson 
with mz' - 1 TeV are 37%, 11%, and 9%. The systematic uncertainties on the signal efficiency vary from 
3% to 6% depending on the resonance mass and decay mode. The primary contributions are due to the 
number of MC events, and the uncertainties related to the muon and t pj scales. 

The observed and expected event yields in both the validation and search mtp regions for all three final 
states are in good agreement, as summarized in Table 1. The mw distributions (Fig. 1) show no significant 
excess above the SM expectation in any of the three modes. The dominant contributions to the uncertainty 
bands in Fig. 1 are due to the number of MC events, the MC cross-section uncertainties, the F™**® scale 
and resolution, and the uncertainty in the shape of the mw distribution for IT-i-jets. 

Upper limits are placed on the production cross section times branching ratio [cr{pp —» v^/Z')xBR(Vt-/Z' ^ 
££'y]. For each or Z' mass, m, the search region is defined to be m + 3cr^r', where o-a> is the standard 
deviation of the simulated signal ma> distribution. The relative width of the signal rriic distribution ranges 
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Table 1: Estimated SM backgrounds and observed event yields for the validation {niw < 200 GeV) and search 
{mw > 200 GeV) regions. Both the statistical and systematic uncertainties are included. Due to correlations, the 
total uncertainties are not exactly the sum in quadrature of the components. 



mu’ < 200 GeV 

ma’ > 200 GeV 

Process 

Net, 

■Verhad 

■V/iThad 

Nef, 

■VeThad 

■Vprhad 

Z/y* —> TT 

6000±400 

11000± 900 

11200± 700 

28± 12 

72+ 21 

99+ 33 

Z/y* —> ee 

— 

6100±1100 

— 

— 

430+ 70 

— 

Z/y* niu 

— 

— 

19500±1300 

— 

— 

410+ 80 

tt 

4220±290 

690± 60 

580± 50 

1640±120 

700+ 60 

550+ 40 

Diboson 

1440± 80 

321± 29 

258± 17 

474± 30 

197+ 17 

141+ 11 

Single top quark 

470± 40 

87± 11 

60± 7 

202± 17 

90+ 10 

73+ 8 

IT4-jets 

54± 18 

17000±4000 

14000±4000 

+1 

CO 

3600+700 

2800+600 

Multijet 

227± 32 

4800±1000 

700± 800 

58+ 12 

340+210 

100+190 

Total 

12400±600 

40400±2900 

46000±4000 

2400+130 

5400+500 

4200+400 

Data 

12954 

41304 

48304 

2474 

5336 

4184 


from 3% to 17% for different mass points, channels and models. To reduce the statistical error, if the up¬ 
per side of the search region is greater than 1 TeV, all events above 1 TeV are used. To further reduce 
the effect of fluctuations in the high-mass region due to low MC event counts, the number of background 
events in each mass window is estimated using a double exponential fit to the total background mw dis¬ 
tribution. The fit uncertainty is taken into account in the limit-setting procedure, including a contribution 
from varying the fit function range. 
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Figure 1; Observed and predicted efi, erhad, jurhad invariant mass distributions. The contributions of the different 
processes are also shown: “Others” includes diboson and single-top while “Jet fake” refers to VT4-jets and multijet. 
All overflows are included in the rightmost bin. Signal simulations are shown for mj,^- 1 TeV and mz> - 0.75 TeV. 
The couplings = 0.11 and Ask - 0.07 (Qu’ = 1) are used for the RPV (Z') model. The uncertainty bands 
include both the statistical and systematic uncertainties. 


A frequentist technique [38] is used to set the expected and observed upper limits as a function of niy^ and 
mz '. The likelihood of observing the number of events in data as a function of the expected number of 
signal and background events is constructed from a Poisson distribution for each ma' bin. Systematic un¬ 
certainties are taken into account with Gaussian-distributed nuisance parameters. A 95% confidence level 
(CL) limit is then determined. The expected exclusion limits are determined, using simulated pseudoex- 
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periments containing only background processes, as the median of the 95% CL limit distributions for each 
set of pseudoexperiments at each value of niy^ or mz’, including systematic uncertainties. The ensemble 
of limits is also used to assess the Icr and 2cr uncertainty envelopes of the expected limits. 

Figure 2 shows the observed and expected cross section times branching ratio limits as a function of 
rriy^imz'), together with the ±lcr and +2cr uncertainty bands. For a Vr mass of 1 TeV, the observed limits 
on the production cross section times branching ratio are 0.5 fb, 2.7 fb, and 9.1 fb for the efi, er and jut 
channels, respectively. The corresponding limits for a Z' boson mass of 1 TeV are 1.0 fb, 4.0 fb and 9.9 
fb for the e/u, er and /it channels, respectively. 

Theoretical predictions of cross section times branching ratio [33,39] are also shown, assuming T'jj = 
0.11 and Ask - 0.07 for the Vt and Qij = 1 for the Z', consistent with benchmark couplings used in 
previous searches. 

For these benchmark couplings, the lower limits on the Vt mass are 2.0 TeV, 1.7 TeV, and 1.7 TeV for the 
efi, er and /it channels, respectively. The corresponding lower limits on the Z' mass are 2.5 TeV, 2.2 TeV 
and 2.2 TeV for the c/i, er and /it channels, respectively. The observed lower mass limits are a factor 
of three to four higher than the best limits from the Tevatron [13, 14] and also more stringent than the 
previous limits from ATLAS [15] for the same couplings. 

In summary, a search has been performed for a heavy particle decaying to e/u, eThad^ or /iThad final states 
using 20.3 fb“' of pp collision data at = 8 TeV recorded by the ATLAS detector at the LHC. The data 
are found to be consistent with SM predictions. Limits are placed on the cross section times branching 
ratio for an RPV SUSY Vy and a LFV Z' boson. These results considerably extend previous constraints 
from the Tevatron and LHC experiments. 



Figure 2: The 95% CL limits on cross section times branching ratio as a function of U mass (top plots) and Z' mass 
(bottom plots) for ep (left), er (middle), and pr (right). Theory curves are for the arbitrary choice of couplings 
djjj = 0.11 and Aj^k - 0.07 for U and - 1 for Z'. The gray band around the theory curve represents the 
theoretical uncertainty from the PDFs and factorization and renormalization scales. 
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